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Abstract 

A sequence inverse relationship can be defined by a pair of infinite inverse matrices. 
If the pair of matrices are the same, they define a dual relationship. Here presented 
is a unified approach to construct dual relationships via pseudo-involution of Riordan 
arrays. Then we give four dual relationships for Bernoulli numbers and Euler num¬ 
bers, from which the corresponding dual sequences of Bernoulli polynomials and Euler 
polynomials are constructed. Some applications in the construction of identities of 
Bernoulli numbers and polynomials and Euler numbers and polynomials are discussed 
based on the dual relationships. 
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1 Introduction 

Krattenthaler defines a class of matrix inverses in [26], which has numerous famous special 
cases presented by Gould and Hsu [13], Carlitz [4], Bressoud [2], Chu and Hsu [9], Ma[29], 
etc. Hsu, Shiue, and one of the authors study Riordan matrix inverses in [22]. Let infinite 
low triangle matrices (d n ,k)o<k<n and ( d n ,k)o<k<n be matrix inverses. Then a sequence 
inverse relationship can be defined as 

n n 

fn = ^ ^ d n ^k3k ' f 7 g n = y ] d n ^kfki (1) 

k—0 k =0 

where sequences {f n }n> o and {g n }n >o are called the inverse sequences with respect to 
inverse matrices ( d n>k )o<k<n and (d nyk )o<k<n- If d n , k = d n , k , he., {d n , k }o<k<n is self¬ 
inverse, then {f n }n> o and {g„}n >o are said to be a pair of dual sequences (or they are dual 
each other) with respect to self-inverse matrix (d n ,fc)o<fc<n, he., they satisfy 


1 
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fn — ^ ^ dn,kQk ^ ^ 9n — ^ ^ ^n,fc/fc? (2) 

/c—0 /c—0 

If there exists a sequence {f n }n >o that is dual to itself with respect to a self-inverse matrix 
(dn,/c)o<fc<ri 5 i.e., 


fn=J2<kfk, 


( 3 ) 


k—0 


then, {f n }n >o is called self-dual sequence with respect to the matrix ( d n ,k)o<k<n • 
Let d nt k = (fc)(-l) fe - Then (d„, fe ) o<fe<n is a self-inverse matrix because 


ee(;)G)(-d‘ + j =^. 

fc=Oj=0 v ' yJ/ 

where S n j is the Kronecker symbol. Let {a n }n>o be a sequence of complex numbers. Then 
«}n>o defined by 


“» = EftV 1)% (<) 

k—0 ' ' 

is the dual sequence of {a„} n >o (see, for example, Graham, Knuth, and Patashnik [14]) 
with respect to ((")(—l) fe )- Hence, 


a n = X (t) ( - 1 ) fca fe- ( 5 ) 

k =0 ' ' 

Wn}n >o an d {a)t}n>o are a pair of inverse sequences with (a*)* = a n . If a* = a n , then 
{a n }n>o is called the self-dual sequence (see Z. W. Sun [48]). For instance, the following 
number sequences are self-dual sequences with respect to the dual relationship (4) (see, for 
example, Z. H. Sun [47]): 

|n + 2m—1| ^ > {( — 1) B n }, {L n }, {nF n - 1}, 

where {B n }, {L n }, and {F n } are the Bernoulli sequence, Lucas sequence, and Fibonacci 
sequence. 

Bernoulli polynomials B n (x) and Euler polynomials E n (x) for n = 0,1,... are defined 
by 



te xt 
e* - 1 


E +n O p xt 

B n (x)- and 
n! 

n>0 


e* + 1 


^2 E n(x) 
n> 0 


t n 

n\ 


Bernoulli numbers B n and Euler numbers E n for n = 0,1,... are defined by 


( 6 ) 


Bn ■= B n ( 0) and E n = 2 n E n 


1 

2 


( 7 ) 
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A large literature scatters widely in books and journals on Bernoulli numbers B n , and 
Bernoulli polynomials B n (x ) and Euler numbers E n and Euler polynomials E n {x). They 
can be studied by means of the binomial expressions connecting them, 


B n (x) = Y,(^)B k x n - k , n> 0, (8) 

^(*) = £(3Mr fc fr. (9) 

fc=o ' ' ' ' 

where E}. = 2 fc E/.(l/2). The study brings consistent attention of researchers working in 
combinatorics, number theory, etc. 

This paper will study the duals of Bernoulli number sequence and Euler number sequence 
with respect to ((]])(—l) fc ) and other self-inverse matrices. The dual sequence of {-B„}n>o 
with respect to ((£)(—l) fc ) is denoted by {B *} n > 0 , and the corresponding dual Bernoulli 
polynomials denoted by o is defined similarly to (8) by using the the duals of 

sequences {B n }n>o, he., 


K(x) = p^ B * k x n ~ k , n > 0 , ( 10 ) 

where R* are the duals of the Bernoulli numbers {B n } n > 0 ! i.e., 

B - = E(fcV 1)fci4 ’ n -°- (11) 

k =0 ' ' 

Hence, there are several questions raised: (1) With respect to which self-inverse matrix 
the Bernoulli number sequence {B n } is self-dual? (2) What is the relationship between the 
self-inverse matrix in (1) and the self-inverse matrix for {(—l) n B n } n >o, ((£)(-l) fc )? And 
(3) does there exist a unified approach to construct self-inverse matrices and self-duals with 
respect to a self-inverse matrix? We will answer those questions by using the Riordan array 
theory. 

Riordan arrays are infinite, lower triangular matrices defined by the generating function 
of their columns. They form a group, called the Riordan group (see Shapiro et al. [44]). 
Some of the main results on the Riordan group and its application to combinatorial sums 
and identities can be found in [5] [8], [10] [12], [15]—[25], [27]-[28], [31]-[32], [34], [36]-[37], 
[39]-[43], [45]—[46], and [50]-[53], 

More formally, let us consider the set of formal power series (f.p.s.) E = R[f]; the order 
of /(f) e E, f(t) = YlkLo fkt k (fk G M), is the minimal number r € N such that f r / 0; T r 
is the set of formal power series of order r. It is known that Eq is the set of invertible f.p.s. 
and E\ is the set of compositionally invertible f.p.s., that is, the f.p.s. fit) for which the 
compositional inverse /(f) exists such that /(/(f)) = /(/(f)) = t. Let d(t) € -Ao and h(t) € 
Ei, the pair (d(t), h(t)) defines the (proper) Riordan array D = (d nj fe)n,feeN = (■ d(t),h(t )) 
having 

d n ,k = [t n ]d(t)h(t) k 


( 12 ) 
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or, in other words, having d(t)h(t) k as the generating function whose coefficients make-up 
the entries of column k. 

It is immediately to be known that the usual row-by-column product of two Riordan 
arrays is also a Riordan array: 

/ii(t)) * (d 2 (t), h 2 (t )) = (d 1 (t.)d 2 {h 1 (t)), (*)))> (13) 

where the fundamental theorem of Riordan arrays 

(d(t), h(t))f(t ) = d(t)f(h(t)) (14) 

is applied to obtain (13). The Riordan array I = (1, t) is everywhere 0 except that it 
contains all l’s on the main diagonal; it can be easily proved that I acts as an identity for 
this product, that is, (1, t) * (d(t), hit)) = ( d(t ), h(t)) * (l,t) = (d(t), h(t)). From these 
facts, we deduce a formula for the inverse Riordan array: 

Mt)) ” = (w i(,) ) (15) 

where h{t) is the compositional inverse of h{t). In this way, the set 1Z of proper Riordan 
arrays forms a group. 

From [40], an infinite lower triangular array [dn,fc]n,fceN = ( d(t),h(t )) is a Riordan array 
if and only if a sequence A = (ao ^ 0, aq, a 2 ,. ..) exists such that for every n, k G N there 
holds 

d n +i,k+i = cio d n ,k + aid n ,k+i + • • • + a n d n<n , (16) 

which is shown in [24] to be equivalent to 

h(t) = tA(h(t)). (17) 

Here, A{t) is the generating function of the H-sequence. In [24, 31] it is also shown that 
a unique sequence Z = {zo, z\, z 2 ,...) exists such that every element in column 0 can be 
expressed as the linear combination 


d n + 1,0 — Zodnfi + Z\dn t l + • • • + Z n d n ^ n , 


(18) 


or eciuivalently, 


d{t) 


dp,o 

1 — tZ{h(t))' 


(19) 


We call sequences A and Z the A- (characterization) sequence and Z— (characterization) 
sequence of the Riordan array ( d{t),h(t )), respectively, and the generating functions of 
called H-sequence and Z —sequence the A (generating) function and Z (generating) function, 
respectively. 

In next section, we present a unified approach to construct a class of self-inverse matrices 
and their application for the construction of dual number sequences and dual polynomial 
sequences. In Section 3 we will discuss the algebraic structure of dual number sequences 
with respect to the dual relationships established in Section 2. Some identities of self-dual 
number sequences are found accordingly. In Section 4, more identities of dual number 
sequences will be constructed by using the dual relationships. 
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2 Construction of self-inverse matrices and the corre¬ 
sponding self-duals 

It is calear that a Riordan array (d(t), h(t)) and its inverse 

(d(t), hit))- 1 = ^ - 1 ^ , h(t)^j , where h(h(t)) = h(h(t)) = t, 

is a pair of inverse matrices. If (d(t), h{t))~* = ( d(t),h(t )), i.e., (d(t),h(t)) is an involution, 

1. e., it has an order of 2, then (d(t), h(t)) is a self-inverse matrix. Here, if g is an element of 
a group G, then the smallest positive integer n such that g n = e, the identity of the group, 
if it exists, is called the order of g. If there is no such integer, then g is said to have infinite 
order. It is well-known that (see Shapiro [42]) if we restrict all entries of a Riordan array 
to be integers, then any element of finite order in the Riordan group must have order 1 or 

2, and each clement of order 2 generates a subgroup of order 2. Sprugnoli and one of the 
author find the characterization of Riordan arrays of order 2 in [24]. 

In combinatorial situations, a Riordan array often has nonnegative integer entries and 
hence it can not have order 2. Therefore, we consider (see definition, for example, in [6]) 
an element R £ 1Z to have pseudo-order 2, i.e., RM has order 2, where M = (1, —t). Those 
R are called pseudo-Riordan involutions or briefly pseudo-involutions (see Cameron and 
Nkwanta [3] and [42]). We now present some sufficient and necessary conditions to identify 
pseudo-involutions. 

Theorem 2.1 Let (d(t),h(t)) be a pseudo-involution, and let A{t) and Z{t) be the gen¬ 
erating functions of the A-sequence and Z-sequence of (d{t),h(t)), i.e., A-function and 
Z—function. Then the following statements are equivalent to the statement that the Ri¬ 
ordan array ( d(t),h(t)) is a pseudo-involution: 

(1) ±(d(t), h(t))(l, — t) = i±d(t), —h(t)) are involutions. 

(2) (l,—t)(d(t),h(t))(l,—t) = (d{t), hit)) -1 , the inverse of (d(t), h{t)). 

(3) ±(1, — t)id(t), hit)) = (±d(— t), h(—t)) are involutions. 

U) = T0i) and z (t) = ^ir- 

Proof. Let (ci(f),h(f)) be a pseudo-involution. Then (d(f), /i(t))(l, — t) = (d(f), — hit)) is an 
involution, i.e., 


id{t),-h{t))id(t),-h{t)) = id{t)di-h{t)),-hi-h(t)) = (1 ,t). 

Additionally, — (d(t), h(t))(l, — t) = (— d)t), — h)t)) satisfies 

i-d(t),-h{t))i-d{t),-hit)) = id{t),-h{t)){d{t),-h{t)) = (l,i), 

which implies — (d(f), h(t))(l, — t) = (— d(t), — hit)) is also an involution. Hence, we have 
proved the assumption that (d(f), h(t)) is a pseudo-involution implies (1). (2) follows from 
(1) because 


/ = Udit), h(t))il, -t)) Udit), h(t))il,-t)) = idit), hit.)) ((1, -t)(d(t), h(t))(l, -t)). 


From the above equations, we have 
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(±(1, -t){d(t),h{t))) (±(1, -t)(d(t),h(t))) 

= ((!, h(t))( 1, -t)) ( d{t),h(t )) = (1, t) = I. 

Hence, we obtain (3) from (2). Similarly, (1) follows from (3). 

To find the A-function and Z-function of pseudo-involution (d(t),h(t)), we recall (see 
Theorem 3.3 of [24]) the A-sequence of 

(d 3 (t),h 3 (t)) = (di(t), h!(t))(d 2 (t), h 2 (t)) 

has the A- (generating) function 

A 3 (t) = A 2 {t)A x ) ’ 

where A\(t) and A 2 (t) are the A- (generating) functions of (di(t), hi(t)) and (d 2 (t), h 2 (t)), 
respectively. Since the H-functions of ( d(t),h(t )) and (1, —£) are respectively A(t) and —1, 
the A-function of (d(t),h(t))(l,—t) is 


A 3 (t) = -A(-t). (20) 

On the other hand, from Theorem 4.3 of [24] we know the A-(generating) function of the 
involution (d(t),h(t))(l,—t) = ( d(t),—h(t )) is 


A 3 (<) = 


-h(t )' 


( 21 ) 


Comparing (20) and (21) we have 


m = 


—t 

h(—t )' 


From Theorem 3.4 of [24], the Z- (generating) function of 

{d 3 {t),h 3 (t)) = ( d(t ), h(t))( 1, -t) = (d(t), —h(t)) 


IS 


Z 3 (t) = Z(—t) 

because the Z-functions of ( d(t),h(t )) and (1 ,—t) are Z{t) and 0, respectively. Since 
(d 3 (t), h 3 (t)) is an involution, from Theorem 4.3 of [24] we also have 


z 3 (t) 


1 — d(t) 
- h(t ) 


Comparing the last two equations, we immediately know 


Z(t) 


d(—t) — 1 
h.(-t) 


From the definition of A-sequence and Z-sequence, a Riordan array (d(t), h(t )) that possesses 
the above A-function and Z-function is a pseudo-involution. 
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Corollary 2.2 Let nonzero d(t) £ T o and h(t) £ T\. Then the infinite lower triangle 
matrix ( d(t ), h(t)) is a pseudo-involution if and only if 

h{t) = -h(-t) and At)= h( _ t) h }~fl_ t)+t , 
where the denominator is assumed not being zero. 

Proof. h(t) = —h(—t) if and only if 


t —t 
h(t) h(—t) ’ 

or equivalently, there exists the function Aft) such that 

A{t) = W) and #) = ^t)’ 

which implies that (d(t),h(t)) has an A- function (see, for example, [24]) satisfying (4) of 
Theorem 2.1. 

d(t) = h(—t)/(h(—t) — td(—t) +t) if and only if 

d(—t) — 1 _ d(t) — 1 
h(—t) td(t) 

or equivalently, there exists the function Z (t) such that 


Z(t) 


d(t) — 1 
td(t) 


and Z{t) 


d(—t ) — 1 
h(~t ) 


which implies that ( d(t),h(t )) has a ^-function (see, for example, [24]) which satisfies (4) 
of Theorem 2.1. Combining all above statements together we know (d(t), h(t)) is a pseudo¬ 
involution, completing the proof of Corollary 2.2. 


Corollary 2.2 also presents an algorithm to find a pseudo-involution. Generally one can 
accomplish this by carrying through the procedure demonstrated by the following example 
as suggested by the corollary. For instance, it is clear that h(t) = t/( 1 — t) (or t/{ 1 + 1)) 
has the compositional inverse h(t) = t/( 1 + t) (or t/( 1 — t)) and satisfies h(t) = —h(—t ). 
From Corollary 2.2, we may also find that d(t) = 1/(1 — t) (or 1/(1 + t)) satisfies d(t) = 
h(—t)/(h(—t) — td(—t) + t). Therefore (1/(1 — t),t/( 1 — t)) ((1/(1 + t),t/( 1 + t))) is a 
pseudo-involution. 

Remark It can be seen that ( d(t),h(t )) is a pseudo-involution if and only if there exist 
^4-function and ^-function satisfying four relationships 


Am = Wr m = wry z <*> 


d(t) — 1 
td(t ) 


, and Z{t) 


d(—t) — 1 
h{—t) 


which not only give the formulas shown in Corollary 2.3 but also the following useful formula 
in the construction of pseudo-involutions: 
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h(t) 


tZ(t) 

Z(-t)(l-tZ{t)) 


and d(t) 


h(t)Z(—t ) 
tZ(t) 


The proof is straightforward from the above four relationships and is omitted. 
From (1) and (3) of Theorem 2.1 we obtain the following result. 


Corollary 2.3 {1/{l-t), t/{l-t)) = {{l)) n k > 0 and {l/{l+t),t/{l+t)) = (©(-1)" fe )„, fc >o 
are pseudo-involutions. Hence, from (1) and (3) of Theorem 2.1, we generate the following 
four Riordan involutions, denoted by R±, R 2 , R 3 , and R 4 , respectively. 


Ri 


i?2 


Rs 


R 4 


1 t 


l-t l-t 
1 -t 


l-t l-t 
1 t 




(- 1 )' 


1 t 


1 + f’ 1 + t 


n,k> 0 


i-t’ l-t 

1 t 


(i,-t) = -(i,-t) 


1 t 


i + t’ i + t 


t - 1 t - 1 

1 t 


(-1) 


fc+i 


(i, -*) 


i-t i-t 

1 -t 


n,k> 0 

1 t 


i + t 1 1 + t 


(i ,-*) 


i+t i+t 


(-1 r 


-(i ,-*) 


1 t 


i-t i-t 
1 -t 


n,k> 0 

1 t 


1 + t’ i + t 


(i ,-t) 


1 H - t 1 -\-1 


(-1) 


n+1 


n,k> 0 


( 22 ) 


Theorem 2.4 iet (i = 1,2,3,4^1 be the Riordan arrays shown in Corollary 2.3. Then 
there hold the following four dual relationships, denoted by D 1 , D 2 , D 3 , and D 4 , respectively. 


Di 


: a n 





r >2 : a n 



(-1 ) k+ 1 a k , 


D 3 : a n 



(-1 ) n a k , 


D 4 : a n 



(-1 ) n+ 1 a k . 
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Furthermore, {(— l) n B n } n >o and {-B ra }n>o are self-dual sequences with respect to D\ and 
D%, respectively. { E n (1) — and {(— l) n ( E n (^) — ^r)} are self-dual sequences with 
respect to D 2 and D 4 , respectively. 

Proof. Since (see, for example, Apostol [1] and Milton and Stegun [33]) 

Y j ( n \B k = B n {l) = {-l) n B n , (23) 

k —0 ' ' 

{(— l) n B n } n > 0 and {B n } n >0 are self-dual sequences with respect to D\ and D 3 , respectively. 
Similarly, from (9) and (—l) n E n (—x) = —E n (x) + 2x n (see, for example [1, 33]) there hold 


k—0 


E I (-D 


\fc+i 


Ek 




fc =0 


1 1 

'2 ~~ 2 


1 

” 2 f 

n—k 


Ek 


= E n -2— + — =E n - — 

\ 2 1 2 n 2 n I 2/ 2” 


k—0 


n \ / 1 


^U]*' ^ 2 fc 


Similarly, we have 


E(:)(-ir'(-i)‘(^a)-F 


k =0 




k—0 

~ E " (-5) + (- 1>n w 

(- 1 )" ( £ " G) - f) ■ 


fc =0 


2 fe 


completing the proof of the theorem. 


We now consider the duals of Bernoulli and Euler numbers and the corresponding duals 
of Bernoulli and Euler polynomials with respect to different dual relationships shown in 
Theorem 2.4. 

Theorem 2.5 Let be the duals of Bernoulli numbers B n with respect to R \, and let 
Bif{x) be the corresponding dual Bernoulli polynomials defined by (10). Then there hold 


B* n (x) = (-l) n B n (-x-l) 


( 24 ) 
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for all n > 0, and 


Proof. (8) and (10) give 


B* n = (-irB n +n. 


Hence, 




k—O 


= e : (*r- fe £ : (-d^- 


!=o 


j=o 


k=j 


- E(")(-o^E(n)w 


k-j 


i—k 


E(”) ( - i) ^i ( ^+ i r- j 

1=0 

(—1)™J5„(—a: — 1). 


which implies (25). 


B; = s;(o) = (-i) n B n (-i) 

B n (l)+n = (-l) n B n (0)+n, 


(25) 


Corollary 2.6 Lei B^(x) be the duals of Bernoulli polynomials. Then their generating 
function is 


j-7i j.n 

n>0 n>0 

_^ e (^+ 1 )i e (x+l)t 

e-* - 1 _ i + ’ 

Other duals of Bernoulli and Euler numbers can be presented below. 


(26) 


Theorem 2.7 With respect to the dual relationship D 3 , the duals of numbers (—1 ) n B n , 
denoted by ((—1 ) n B n )*, may be written as 

((-1 YBnY = £ (—1)"( — l) fc -Bfe = B n + (—1 )"n. 

And the corresponding dual Bernoulli polynomials are 


(27) 
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E (?) ((-1 ) k B k )*x n ~ k = B n (x) - n(x - I)"" 1 . (28) 

k—O ' ' 

With respect to the dual relationship D 4 , the duals of numbers E n (l/2 ) — (l/2) n , denoted 
by (E n ( 1/2) — (l/2) n ) , have the expressions 



The corresponding dual Euler polynomials are 


E Q (*- k Q) - = {-l) n E n {-x+l)-2{x-l) n + {x-2y. (30) 

With respect to the dual relationship Di, the duals of numbers (—l) n (E n ( 1/2) — (l/2) n ), 
denoted by ((—I)" (E n ( 1/2) — (l/2) n ))*, can be evaluated as 



And the corresponding dual Bernoulli polynomials are 

± 0 (*-(<-.)* © - ^))* -ft W +<. + 1)- - *•. <») 

Proof. The proofs are straightforward from the definitions and are omitted. 


3 Generating functions of self-dual number sequences 

In this section, we give some structures of self-dual number sequences, by using which the 
characterizations, relationships, and other properties of the self-dual number sequences can 
be obtained. The first half of the following theorem is presented in Prodinger [38]. 

Theorem 3.1 Let a(x) = Z] n >o a nX n - Then its coefficient sequence is a self-dual sequence 
with respect to D\ (D 2 ), i.e., it satisfies 
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for ngNU {0} if and only if a(x) satisfies the equation 


a = (! - x)a{x) (-(1 - x)a{x)), 

respectively. The coefficient sequence of a{x) is a self-dual sequence with respect to (D 4 ), 

i.e., it satisfies 


^ (-1 ) n a k = a n (—a n ) 

k =0 ' ' 

for ngNU {0} if and only if a(x) satisfies the equation 

= (l+x)a(x) (— (1 + x)a(x)), 

respectively. 

Proof. Since the first half of the theorem is given in [38]. It is sufficient to prove the second 
half. 



±- 


1 + x V 1 + x 


= ±^Ta fc (-l)V=(l + : 


\—k— 1 


= ±^a fe (-l) fe x fc ^ 
k =0 j =0 


k —0 

-k- 1 
j 


OO OO 


k =0 j =0 

00 / n / 

= ±e 


k + j 
j 


k+j 


n—0 


3=0 


an- 7 X 


= ± a n x n = ±a(x). 


n—0 


The proof is complete. 


Corollary 3.2 Let {a n } n >0 be a given sequence with ordinary generating function a(x) = 
12 n >o a nX n . Then 

(1) [47] {a„}„>o is a self-dual sequence with respect to D\ if and only if {2a n +i — a n } n >o 
is a self-dual sequence with respect to D 2 . 

( 2 ) {a„} n >o is a self-dual sequence with respect to D 3 if and only if {2a n+ i + a n } n >o is 
a self-dual sequence with respect to D 4 . 
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Proof. The result shown in (1) is given in [47]. It is sufficient to prove (2). Let 

b n = 2a n +i + a n 

with its ordinary generating function b(x) = bnX n . Thus, 

2(a(x) -o 0 ) , , ^ x + 2 2 

blx) = -b a(x) = - a(x) - a o, 

x x x 

which implies 


Therefore, 


if and only if 


b -- 


1 + x 


x + 2 / x 

- a -- 


, i x + 2 1 

-(1 +x) | - ——a | - 


X 1 + X V 1 + x 


1 + X 
X 


H — (1 H - x)clq 

X 


- a 0 • 

X 


1 + X 


= (1 + x)a(x) 


b ( - X ^ = -(1 + x) ( X + 2 a(x) - — a 0 ) = -(1 + x)6(x). 
1 + x / V ^ x 


Based on Theorem 3.1 we have finished the proof. 


Theorem 3.3 Let {a n } n >o be a given sequence with exponential generating function a*{x) = 
E„> 0 a nX n /n\. Then 

(1) [47] K} „>o is a self-dual sequence with respect to Di if and only if a*[x)e is 
an even function. 

(2) [47] {a n } „>o is a self-dual sequence with respect to D 2 if and only if a*(x)e “Z 2 is 
an odd function. 

(3) K} n >o is a self-dual sequence with respect to D 3 if and only if a*(x)e x ^ 2 is an even 
function. 

(4) {«4 n >0 is a self-dual sequence with respect to D 4 if and only if a* (x)e x ^ 2 is an odd 
function. 


Proof. The results shown in (1) and (2) are given in [47]. We only need to prove (3) and 
(4) here by using a similar argument: 
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fc >0 


j> 0 


3- 


k >0 i>o 

= ElB-D' 


r fc+i 

/c!j! 


n>0 \fc=0 


\ X 

a k —r- 
/ n! 


Hence, if 


then 





a*(—x)e x = a(x) and — a(x), 
respectively, which can be written briefly as 

a*(-x)e ~ x/2 = ±a(x)e x/2 . 

If the case of positive sign on the right-hand side holds, i.e., {«„} is a self-dual sequence 
with respect to D 3 , then a*(x)e x / 2 is an even function; while the negative sign holds, or 
equivalently, {a„} is a self-dual sequence with respect to £>4, then the function a*(x)e x / 2 is 
odd. It is easy to see the sufficiencies of (3) and (4) are also true. This concludes the proof 
of the theorem. 


[47] uses Theorem 3.3 to derive a numerous identities. [30] uses umbral calculus to 
reprove and extend some of them. We now survey their results and extend them to other 
self-dual sequences. 

Theorem 3.4 For any function /, we have 

(1) [4V ELo (fc) (m - (-1 r~ k Ej=o ( *n-k = 0 for n € N U {0} */ {a n } n > 0 

is a self-dual sequence with respect to D\. 

( 2 ) [30] ELo (fc) (/(*) + (-l) n_fc Ej=o a n-k = 0 forn e N U {0} if {a n } n > 0 

is a self-dual sequence with respect to Di. 

( 3 ) ELo (k) (/( fc ) - Ej=o(“ 1 ) n_J (i)/(i)) a n-k = 0/orneNU {0} */ {a n }„>o is a 
self-dual sequence with respect to D 3 . 

(4) ELo (fc) (/( fc ) + E 4 =o(“ 1 ) n_J (i)/(i)) a «-fc = 0 forn e NU {0} */ {a n }„>o is a 

self-dual sequence with respect to D 4 . 

Proof. The proofs of (3) and (4) are similar as (1) and (2) by using either Theorem 3.3 or 
umbral calculus. Hence, we omit them. 
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From Theorem 2.4, {(— l) n B n } n >o and {B n } n >o are self-dual sequences with respect to 
D i and D 3 , respectively. {E n (|) — ^7} and {(—l) n (E n (|) — 7^) } are self-dual sequences 
with respect to D 2 and £> 4 , respectively. Therefore, we may use Theorem 3.4 to obtain the 
following identities. 

Theorem 3.5 Let B n and E n be Bernoulli numbers and Euler numbers, respectively. For 
any function f, we have 



E(:)(/(‘)-E(-)-G)/0'))b„.^„, 

£© ((-ir-vw+B-D-Q/o-)) (b^Q) - 5 ^) =0. ,33, 
where the first identity is given in [47]. 

4 Applications of dual sequences to Bernoulli and Euler 
polynomials 

Let a € R. Denote 


C n ,a{x) = (~l) n X! ( X ~ k and G n,a( x ) 


r)a* k (x-a) n 


Then we have the following identities about C nrCt (x) and C* a {x), which is an extension of 
the main results shown in Sun [49]. 


Theorem 4.1 Let C n , a {x) and C* a (x) be defined as (34), and let k,£ £ NU {0} and 
x + y + z = 1 + 2a. Then there hold 



^ 2 (-l) j x e ~ j 


j —0 



G k+j + l,a( Z ) 

k+j + 1 


( 35 ) 
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In addition, we have 


j2(-iyx k -j 


3=0 


C e+j , a (y) = J 2 (- i y x *~ j 

3=0 


C* k+j , a (z) 


and 


E(-i vv+j+ l > k ~ j+1 ( k+c i+j Av) ++ 3 + v 

3=0 \ 3 / j — 0 

= (k + £ + 2) ((— l) k C'k + e + i ta (y) + (—1 Ycf. +e+la (z)). 

Proof. Substituting (34) into the left-hand side of (35) yields 


" fk\ (-1 

~" 0 \jJ V + J + 1 


(v - «)' +>+i + e »,« rtit 1 ) 


i =0 


= ca a + E^EC =Vt + } + 1V ~ 1>V ~’ 


i =0 j—O 

fc+^+1 £ k+j+1 

+ E a 'E E 

i—1 j—O r=i 
k+£ 

= m o + E— 7-ra i+1 , 

l + 1 


jj \ i +1 y^ + j' + i 
A fk+j + 1\ fr\ (—l) j x e ~ j 


ij k + j + 1 


(-1 y +3+1 (y - a) e+j+1 - 

( - ( Z - a)) k+j+1 - 


i =0 


where 


c = 



£ + j + l 


(-1 ) e+ i +1 (y - a) e+j+1 



k+j + 1 


(-(* 


- a)) fe+i+1 " 


(36) 

--f ;>%•>) 

(37) 

(- 1)^ +1 (y-a) /+J ‘- 

- a) fc+j+1 “ r 
-(*+i) 


which will be evaluated later, and 
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1—r 




= 0 , 


where the last step is due to the expression in the last parenthesis is zero following from 
Lemma 3.1 of [49]. 

Finally, we evaluate c as follows. 




3=0 




U\ rv-ct L / f\ rx+y-ot 

*—,■ I ^ - \ . / t k+j dt 


= <-i) ,+ ‘E ‘ ,+J < i ‘+E( Jt- 1 ) 

j=0 " ,0 j=0 

(■y-a rx+y—a 


ry-a px-ty-a 

= (-i y +1 t\t + x ) k dt + (- i) e t k (t-xYdt 

Jo Jo 

/ px-\-y—a px+y—a \ 

= (— 1 y y J (s — x) e s k ds + J (t — xYt k dt) 

= (— 1 Y [ (t — xYt k dt = (— 1 Y [ (tx — xY(tx) k xdt 
Jo Jo 

[ X (l _ i;Yt k dl - x fc+£+i r ( fc+ 1 ) r (^+ !) 

Jo [ j r(k + £+2) 


_ x k+e+1 




(k+e+ i)( fc +y 
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By taking partial derivative with respect to y on the both sides of (35) and noting 
z = 1 + 2 a — x — y, one may obtain (36). Replacing k and £ by k + 1 and l + 1 in (36) taking 
the partial derivative with respect to y , one may obtain (37). The proof is complete. 


The following corollary of Theorem 4.1 is equivalent to the results shown in Theorems 
1.1 and 1.2 of [49]. 

Corollary 4.2 ([49]) Let C n< o(x) and C* 0 (x) be defined as (34), and let k,£ € NU{0} and 
x + y + z = 1. Then there holds 


j~k-j f k \ c t+j+ifi(y) x f Ll yV-j (?\ C fc+j+i,o( g ) _ a 0 x k+e+1 


Et -1 )'* 

1=0 


l=o 


jj £ + j +1 ^ ~ \j) k+j + 1 (k +£+l)( k + t ) 


(38) 


In addition, we have 
k 


E(-l)^ fc - j (^Q+l,o(2/) = E(- 1 ) j ^“ i (-)^+l,oW 

1=0 l=o 


(39) 


and 


E(-i y^ + j + ( k + l )ce+M + + 3 + iy- j+1 + 

d —n \ 3 / a —n \ 3 / 


3=0 \ ^ / j = 0 

= (k + e + 2 ) ((-i) k c k+e+ 1 ,o(y) + (-1 Yc* k+Uh 0 (z)). 

If a k = Bk or (—l) k Bk, then we have 


k £. 

(_ 1)^+1 Y x k ~ j E±i±i ^ + (_i)fc+! Y x e ~ j ^ B k+j+l(z) 


1=0 
r k-\~e ~(-i 


1=0 


j) k+j + 1 


(k + £+ l)( fc f) 
In addition, we have 

k 

„k-j 


1=0 


(-I)^E^L- ^ + i(y) = (-i) fc E 


1=0 


.)Bfc+iW 

3 z 


(40) 


(41) 


and 


(-i fYY+j+ox k ~ j 

1=0 


l 


Be+j(y) + (—1) > (k + j + l)x 


1-3 


l=o 


— (k + £ + 2 ) ((—l) fe _Bfc_j_f + i(y) + (—l) £ Rfc + ^_(_i( 2 ;)) . 


+ 

j 


B k+j{z ) 
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The conjugate Bernoulli polynomials B n (x) is introduced in [23] by their generating 
function as follows. 


t n 


1 + 


l +t 


-e« ~ |> 


t n =0 

where the first few terms of the conjugate polynomial sequence {B n } n >o are 


(42) 


B 0 (x) = 1, 

B ± (x) = x+ 

B 2 {x) = x 2 + x+ 

0 

S 3 (a;) = a; 3 + -a; 2 + -x + 2, 

191 

B±(x) = a; 4 + 2a; 3 + 5a; 2 + 8a; + 

- . , * 5 4 25 , 2 191 76 

B 5 (x) = x 5 + —x 4 + —ar + 20a; 2 H— —x + etc. 

Zi o Do 

Let {B n (x)} n > o be the conjugate Bernoulli polynomial sequence defined by (42), and let 
the conjugate Bernoulli numbers {B n } n > o be defined by B n = B n (0). Then we may find 
that 


B n (x) = J2( n \ n ~ k B k (43) 

fc=0 ' ' 

for all n > 0. 

We define dual sequence of the conjugate Bernoulli number sequence with respect to 
inverse matrix R\. Then the corresponding dual polynomial sequence is 

K(x):=p(fyx n ~ k B* k (44) 

for all n > 0. From Theorem 4.1, we obtain 

Corollary 4.3 Let k, l £ N U {0} and x + y + z = 1. Then there holds 


( k )x k ~ j (-i y +1 

a=o ' J ' 


Be+j+i{-y) 
i + j + 1 



^fc+j+i( z ) 

k + j + 1 


(k+i+i^ty 


( 45 ) 


where B n [t) are defined by (42), and B^(x) are defined by (44). Also 
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(-i) fc E ( k )(-x) k ~ j B (+j (-y) = (-l)'E ( ]Y-xY~ j Bl +j {-z) (46) 

3 =o v3 ' j=o 


and 


(-i) e J 2 ( k+1 ) i-x ) k ~ J+1 + m +J (- y ) 

+(-i) fc Eft 1 ) + j + (-*) 

3=0 \ 3 3 

= (k + e + 2)((-l) e + 1 B k+e+1 (-y) + (-l) k + 1 B* k+e+1 (-z)). (47) 

In (34), substituting a = 1/2 and 

”‘ = e ‘G)-f and oJ = (- 1 )* ( E * ( 5 ) + • OS) 

where the duals a* k are derived by using (29) of Theorem 2.7, we obtain 


A n {x) = (-1 ) n E n (x) - (- l) n x n and A* n (x) = E n (-x + 1) + (2 - x) n - 2(1 - *)". 
Hence, Theorem 4.1 implies 

Corollary 4.4 Let C n ^ a {x) and C* a (x) be defined as (34) with a = 1/2, and let k,£ G 
N U {0} and x + y + z = 2. For a k and a k shown in (48), there hold 


^(-1 ) e+1 x k ~ j — e .±i±±M. + ^(-1 y x e ~ j E k+j+i(~ z + !) 


3=0 


3=0 


jj k+j + 1 

x+y—l 


ry rx+y rx+y-1 

= (—l) f+1 / t e (t + x) k dt— / t k (x — t) e dt + 2 / t k {x — tydt. 

Jo Jo Jo 


(49) 


In addition, we have 


(-1 Y E x k ~ j ( k ) Ee +j (y) (-l)V(a + y) k 

3=0 

= ^(-l)V^(% fe+i (l-^) + ( a: + 2/ ) fc (-^-2(3 : + y-l) fc (l-^. (50) 

3=0 
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The above identities of polynomial sequences can be used to establish identities of number 
sequences. For instance, if x = 1, y = 0, and 2=1, then (49) yields the number sequence 
identity 




£(-F 

l=o 


A E, 


k+j+l 


(0) 


yjj k + j + l 
If x = 1 and y = z = 1/2, from (49) there holds 


(k + £+ l)( k f) 


g ( -i)'+v-» Q +b-dv- Q 


-A+i+i (1/2) 

k + j + 1 


1 


- (fc + £ + l ) (^y +2S Qa + M+iJ - 2B{^,k + lJ+l 

where B(a,a,b) = f^t a ~ 1 ( 1 — t) b ~ 1 dt is an incomplete beta function. From the last two 
identities, we have 


5>irv-* 
1=0 


e+i^k-j (A Ee+j +i (1/2) - Ee +j+ i(0) 
j) l + j + 1 


Ei-ov-qj 

l=o v/ 


i-j (A Ek+j +1 (1/2) - Ek+j +1 (0) 
j) k+j + l 


— 2 ( S ( -,k + 1,1 + 1 j — B (-,k + 1,£ + 1 


Let A and B be any m x m and n x n square matrices, respectively. Clieon and Kim [6] 
use the notation © for the direct sum of matrices A and B: 


A 0 
0 B 


A®B = 

Hence, we may obtain a matrix form of Lemma 3.2 of Pan and Sun [35] 
Theorem 4.5 Let n be any positive integer, and let B(t) and P[t] be defined as 


(51) 


B{f) := (H 0 (AA(A- 


P[t] := 


in—k 


n,k> 0 


and 


respectively. Denote ( x ) = (1, x, x 2 ,.. .) T , D = diag( 0,1,1/2,...), and [X] 
i.e., 


[ 0 ]© 


J l<k<n 
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which implies 


\X]B(x + y) = P[x]DB(y) + [X](x), 


+ = y M + HnX „ 


l l 


for all n > 0 shown in [35]. 


Proof. Noting that B 0 (y) = 1, the left-hand side of (52) can be written as 

LHS of (52) 

= [X]P[x]B{y) 

' 1 

x+(\)B! (y) 

x 2 + (l)xB 1 (y) + (2 )B 2 (y) 

= [*] : 

X n + (n) x n-lB 1 (y) + ... + ^)B n (y) 



■ 0 

(DxB^ + QB^y) 
^) x ^B 1 (y) + -.- + OB n (y) 


(1 + ^ 2 


+ [X]\ : 


(^)xBi(y) + (l)B 2 (y) 


(! +1 + • • • + (?)®"- 1 B 1 (y) + • • • + (C)B4y) 
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where the second term can be written as 

r o 


(§+\§)Bi{y)x + \§B 2 (y) 

((!) + *(?) + ■■• + *(?)) + (3(2)+ 5(2) +••• + ^(2)) W”" 2 + • • • + 

' 0 

(l)Bi(y)x+ \ ( 2 ) B 2 (y) 

(^(y)®"" 1 + i ( n 2 )B 2 (y)x^ + • • ^0^(2/) 


In the last step we use the identities i(?) = and 


n - 1\ In - 1\ In - 1\ In - 2\ /n - 3 
ft-l) + ( ft )-(fc-l) + (fc-l) + ( ft 


n — 1\ In —2 

ft - 1/ + Ift- 1 


ft 

ft -1 


for ft,n € N with 1 < ft < n. Since the last matrix can be written as P[x]DB(y), we obtain 
(52). By multiplying the matrices of (52), 


LHS 0 / (53) = 11 n x" + ]T (?) = Z (?) + Hn* n , 


1=1 


1=1 


l l 


which is the RHS of (53). 
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